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I. INTRODUCTION 

The Gamow- Teller (GT) nuclear response is a very fer- 
tile source of information about important issues related 
not only to nuclear physics [![, but also to astrophysics 
and particle physics jH, 4] . In the case of unstable nu- 
clei this information is mainly extracted from /? decays, 
where there is a severe restriction due to the Q- ener gy 
limitation. In the case of stable or close to stability nu- 
clei, the GT strength is obtained from charge-exchange 
reactions at intermediate incident energies and forward 
angles 5]. Under these conditions the nuclear states are 
probed at small momentum transfer and the cross section 
becomes proportional to the GT matrix element without 
the energy limitations that characterize /? decays. 

The spin-isospin nuclear properties in 76 Ge and 76 Se 
are among the most extensively studied both theoret- 
ically [fl-QjJ and experimentally [l3l - [l8j . This is due 
to their significance as double- (3 decay partners and the 
implications of this process in the determination of the 
neutrino nature and its absolute mass 0. The study of 
these nuclei is indeed a part of a large experimental pro- 
gram being pursued in the last several years and aimed 
to explore the GT properties at low excitation energies 
of double-/? decay partners by high resolution charge- 
exchange reactions [J], [H, [I9I - I21I ] . 

The present work is motivated by the recent high reso- 
lution charge-exchange experiment 76 Ge( 3 He, t) 76 As [14[ 
that has allowed the unveilling of some remarkable fea- 
tures of this nucleus, which previous charge- exchange ex- 
periments |13[ at much lower resolution were unable to 
identify. In particular the authors of Ref. fl4| reported an 
unusually strong fragmentation of the GT strength that 
was interpreted in terms of possible effects of deforma- 
tion. In Ref. [HI it was also noted a lack of correlation 
among the GT transition strengths feeding the same lev- 
els in 76 As from the two different directions, GT - mea- 
sured in the 76 Ge( 3 He, t) 7e As reaction and GT+ mea- 
sured in the 76 Se(<i, 2 He) 76 As reaction 16j. In view of 
this new experimental information that has become avail- 
able, it is worth reconsidering the theoretical description 



of these nuclei and the role that deformation might play 
to understand the observed features. 

In this work we explore the ability of the deformed 
proton-neutron quasiparticle random phase approxima- 
tion (QRPA) approach to describe together all of this 
rich information available at present that includes i) the 
global properties of the GT response, such as the total 
GT strength as well as the location and strength of the 
GT resonance, ii) the GT strength distribution in the 
low-lying excitation region that contains much more ac- 
curate information, and iii) the two-neutrino double-/? 
(2i//3/3) decay matrix element and the implications of the 
single (3 branches in the 2vj3(3 process. 

The QRPA is one of the most reliable and broadly 
used microscopic approximations for calculating the cor- 
related wave functions involved in (3 and double f3 [22[ 
decay processes, especially after the inclusion of particle- 
hole (ph) and particle-particle (pp) residual interactions. 
The method was first studied in Ref. [23[ to describe the 
(3 strength in spherical nuclei. Subsequent extensions of 
the QRPA method to deal with deformed nuclei were 
done later in Rcfs. [23 - l29j . Deformation effects were 
also studied in the double-/? decay process B, H , fill l30l — 
[32|]. In particular, it has been found [t| l30l| that the 
nuclear matrix elements for the 2v(3f3 process are sup- 
pressed with respect to the spherical case with a reduc- 
tion factor that scales roughly with the deformation dif- 
ference between parent and daughter. This suppression 
mechanism, which is ignored in spherical treatments may 
play an important role in approaching the theoretical es- 
timates to the experiment. 

As it shall be described later, our theoretical approach 
[28|,[29[ is based on a deformed QRPA formalism on top of 
a selfconsistent deformed Hartree-Fock (HF) mean field 
with Skyrme forces and pairing correlations in the BCS 
approximation. In particular, we shall study the depen- 
dence on deformation of the single /? branches that build 
up the double-/? process. 

The paper is organized as follows. In Sec. II, we 
present a brief summary of the theoretical approach used 
to describe the GT properties. Section III contains the 
results obtained for the GT strength distributions, as well 
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as the results for the 2vj3f3 decay, stressing the effect of 
deformation. The summary and conclusions are given in 
Sec. IV. 



II. THEORETICAL APPROACH 

We describe here briefly the theoretical formalism used 
in this work, whose details can be found in Ref. [28|, |29j . 
First, we carry out a selfconsistent deformed HF calcula- 
tion with the effective nucleon-nucleon density-dependent 
Skyrme interaction SLy4 (33| . assuming axial deforma- 
tion and time reversal symmetry [34j . The single-particle 
wave functions are expanded in terms of the eigenstates 
of an axially symmetric harmonic oscillator in cylindri- 
cal coordinates using twelve major shells. Pairing cor- 
relations between like nucleons are included in the BCS 
approximation taking fixed pairing gap parameters for 
protons and neutrons, which are determined phenomeno- 
logically from the odd-even mass differences of neighbor- 
ing nuclei through a symmetric five-term formula involv- 
ing experimental binding energies [35| . The occupation 
probabilities vf of the single particle levels are computed 
at the end of each HF iteration and are then used to cal- 
culate the one-body density and mean field of the next 
iteration, so that one gets new single-particle wave func- 
tions, energies and occupation numbers at each iteration. 
Therefore, the selfconsistent determination of the bind- 
ing energy and deformation includes pairing correlations 
from the beginning. After convergence, the QRPA equa- 
tions are solved on the deformed ground state basis for 
76 Ge and 76 Se to get their GT strength distributions and 
to compute the 2vf3f3 decay matrix element. 

To describe GT excitations in the QRPA we add to the 
quasiparticle mean field a separable spin-isospin residual 
interaction in ph and pp channels. The advantage of us- 
ing separable forces is that the QRPA energy eigenvalue 
problem is reduced to finding the roots of an algebraic 
equation. The ph part is responsible for the position and 
structure of the GT resonance [ID, [H, [28[ . Its coupling 
constant Xph could in principle be obtained in a consis- 
tent way from the same Skyrme energy density functional 
as the HF mean field through the second derivatives of 
the energy functional with respect to the densities and 
averaging the contact interaction over the nuclear vol- 
ume, as it was done in Ref. 28] to study exotic nuclei. 
The pp part consists of a proton-neutron pai ring force 
and it is also introduced as a separable force [26| , The 
coupling constant is usually fitted to the half-lives 
phenomenology [26] . Following the above mentioned pro- 
cedure and taking into account the experience accumu- 
lated in this mass region [36] . we have chosen in this 

work the values Xph = °- 15 MeV and K pp = °- 03 MeV - 
In addition, we will also show the sensitivity of the GT 
strength distributions and Ivfifi nuclear matrix elements 



to the value of coupling constant '. 

The technical details to solve the QRPA equations have 
been described in Refs. [26, 28, 29]. Here we only mention 



that, because of the use of separable residual forces, the 
solutions of the QRPA equations are found by solving 
first a dispersion relation, which is an algebraic equation 
of fourth order in the excitation energy lo. Then, for each 
value of the energy, the GT transition amplitudes in the 
intrinsic frame connecting the ground state |0) to one 
phonon states in the daughter nucleus \ujk), are found to 
be 



where 



E 1 



(i) 

(2) 
(3) 



with 



in terms of the occupation amplitudes for neutrons and 
protons iv.tt Wt 7t = 1 — v„ n ) and the matrix elements of 
the spin operator connecting proton and neutron single- 
particle states, as they come out from the HF+BCS cal- 
culation, and are the forward and backward 
amplitudes of the QRPA phonon operator, respectively. 

Once the intrinsic amplitudes in Eq. (JT|) are calcu- 
lated, the Gamow- Teller strength £>(GT) in the labora- 
tory frame for a transition IiKi(0 + 0) — > IfKf(l + K) can 
be obtained as 



BuiGT*) = ^ [<^=o l^o^l 0> 2 J(o;^ =0 - u;) 

+2 (lj k=1 Ict^Io) 2 ^^! -u) ,(5) 

in [g A /4ir] units. To obtain this expression we have used 
the initial and final states in the laboratory frame ex- 
pressed in terms of the intrinsic states using the Bohr 
and Mottelson factorization [37]. Finally, a quenching 
factor gA,eS = 0.7 gA.frcc is included in the calculations 
to take into account in an effective way all the correla- 
tions [H[ that are not properly considered in the present 
approach. 

The role of the residual interactions and BCS cor- 
relations on the GT strengths was already studied in 
Ref. [HI, Hi| . The role of deformation was also studied 
there, where it was shown that the GT strength distri- 
butions corresponding to deformed nuclei are much more 
fragmented than the corresponding to spherical ones, be- 
cause of the broken degeneracy of the spherical shells. It 
was also shown that the crossing of deformed energy lev- 
els, which depends on the magnitude of the quadrupole 
deformation as well as on the oblate or prolate character, 
may lead to sizable differences between the GT strength 
distributions corresponding to different shapes. These 
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features have been exploited to use the /3-decay prop- 
erties as an alternative method to learn about the nu- 
clear deformation in highly unstable isotopes [3!| . It has 
also been shown that deformation is a key ingredient 
to reproduce the occupation probabilities of the relevant 
single particle levels in the valence shells of 76 Ge and 76 Se 
involved in the double-beta decay process that have been 
measured for neutrons [17| and protons f 1 81 ) - 

The nuclear double-/? decay is a rare second order weak 
interaction process that takes place when the transition 
to the intermediate nucleus is energetically forbidden or 
highly retarded. Two decay modes are expected, the two 
neutrino mode, involving the emission of two electrons 
and two neutrinos, and the neutrinoless mode with no 
neutrino leaving the nucleus. Whereas the first type is 
perfectly compatible with the Standard Model, the sec- 
ond one violates lepton number conservation and implies 
the existence of a massive Majorana neutrino. Because 
the nuclear wave functions and the underlying theory for 
treating the neutrinoless and the two-neutrino modes are 
similar, the Gamow- Teller part that drives the 2vf3(3 de- 
cay provides insight for theoretical models that are re- 
quired to reproduce the available experimental informa- 
tion on the 2v(3/3 half-lives. 

The 2vj3P decay is described in the second order 
perturbation of the weak interaction as two successive 
Gamow- Teller transitions via virtual intermediate 1 + 
states. The basic expressions for the 2vf3fi deca y w ithin a 
deformed QRPA formalism can be found in [91.I40}. Here 
we only write the half-life of the 2vf3(3 decay 



T, 



1/2 



M, 



2vf3f3 
GT 



(6) 



in terms of the phase-space integral G 2v ^^ and the nu- 
clear matrix element that contains all the infor- 
mation of the nuclear structure involved in the process, 



M, 



E E 

K=0,±1 nii,rrif 



1 \K {^K^rrif \^K, mi ) 
K +0J K )/ 2 



X(0f\a- K t \LL>K,m f ) (uK,rn t \cr K t |0j) . 



(7) 



In this equation (wj^ ) are ^ ne QRPA excitation en- 
ergies of the intermediate 1 + states \uJx,mi)(\^K,m f )) 
with respect to the initial (final) nucleus. The indices 
nii, rrif label the 1 + states of the intermediate nucleus. 
The overlaps are needed to take into account the non- 
orthogonality of the intermediate states reached from dif- 
ferent initial |0j) and final |0/) ground states. Their ex- 
pressions can be found in Ref. [9j. 

The various measurements reported for the 2i//3/3 de- 
cay in 76 Ge have been analyzed in Ref. [4lj, where a 
recommended value = (1.5 ± 0.1) x 10 21 yr was 

adopted. Using the phase-space factor [9J G 2 "' 3 ' 3 ( 76 Ge) = 
1.4910 -20 yr^ 1 MeV 2 , we get the experimental nuclear 
matrix elements M, 



TABLE I: Experimental and calculated charge root mean 
square radii r c [fm], intrinsic charge quadrupole moments Q p 
[fm 2 ], and quadrupole deformations J3 for 76 Ge and 76 Se. Ex- 
perimental values for r c are from [Q] . The first experimental 
values for Q p are from while the second values are from 
H. 



76 Ge exp. 4.080/4.127 66(21)/164.1(2.5) 0.10/0.26 

SLy4 4.104 93.85 0.14 

76 Se exp. 4.088/4.162 119(25) /205.5(2.4) 0.16/0.31 

SLy4 4.151 125.1 0.17 
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quenched factors qa — 1 are used. 
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GT 



0.129 McV" 1 when the bare 



III. RESULTS AND DISCUSSION 

In this section, we first discuss the global ground-state 
properties of 76 Ge and 76 Se. In Table U we compare with 
experiment our results from microscopic SLy4 HF+BCS 
calculations corresponding to charge root mean square 
radii r c , quadrupole moments Q p , and quadrupole defor- 
mation parameters /?. The experimental values for charge 
radii from Ref. [42} are well reproduced in our calcula- 
tions. The experimental intrinsic quadrupole moments 
quoted correspond to values extracted from Coulomb 
excitation reorientation methods (43j and from electric 
quadrupole transition probabilities B(E2) [44]. Our cal- 
culations produce quadrupole deformations within those 
experimental values. 

Besides the selfconsistent solution that gives us the en- 
ergy minimum and the nuclear shapes at equilibrium, 
we are also interested in the behavior of the energy 
as a function of the deformation. These energy curves 
are obtained by performing constrained HF+BCS cal- 
culations [451 ] . where the HF energy is minimized un- 
der the constraint of keeping the nuclear deformation 
fixed. The energy curves for 76 Ge and 76 Se exhibit two 
local minima at oblate and prolate shapes that are prac- 
tically symmetric with very low energy barriers. The 
profiles of these curves are very shallow and roughly 
have the same energy for quadrupole deformations from 
/3 = —0.2 up to /3 = 0.2. These results are in agreement 
with those from similar calculations using the Gogny- 
D1S finite-range effective interaction [46} . Taking into 
account these characteristics, with oblate, spherical and 
prolate shapes having practically the same energy and 
where small changes in the calculation can lead to differ- 
ent results for the absolute minimum associated to the 
ground state, we have opted to show here results for the 
GT strength distributions at various shapes correspond- 
ing to P = -0.2,0,0.15,0.2,0.3 to study the sensitivity 
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TABLE II: Measured and calculated GT strength in 76 Ge accumulated in various energy regions [MeV]. 



£B(GT-) E0-5 E0-7 E 7 -10 ElO-20 Eo -20 
Thies et al. [14] 1.60(18) / 2.43(32) 

Madey et al. [13] 1.52 4.88 2.58 12.43 19.9 

QRPA (B = 0.15) 4i) 5.5 2.2 10.6 18.3 



TABLE III: Same as in Table HH but for 76 Se. 



E B(GT + ) Ep-2 Eq-5 E5-20 Ep-20 

Grewe et al. [16] 0.54(7) 0.7(2) 0.3(2) 1.0(4) 

Helmer et al. [15] 0.38(3) 0.79(7) 0.64(9) 1.43(16) 
QRPA (/3 = 0.2) 0.15 0.87 0.53 1.4 



of the GT strength distributions to the deformation. 

A. Gamow- Teller strength distributions 

In the upper plot of Fig. Q] we show the measured 
GT strength distribution as a function of the excita- 
tion energy of the daughter nucleus, extracted from 
the charge-exchange reactions 76 Ge(p, n) 76 As [l3| and 
76 Ge( 3 He, t) 76 As [lj]. While the former measurements 
were taken at a rather low resolution, the latter exper- 
iment was performed at a high energy resolution of 30 
keV. The inset shows a more detailed comparison of these 
two measurements in the overlapping energy range below 
5 MeV. The curves correspond to the same distributions 
obtained from a folding procedure using 1 MeV width 
Breit-Wigner functions, so that the original discrete spec- 
trum is transformed into a continuous curve. 

In the lower panel we can see the GT strength distri- 
butions calculated within QRPA for various deformation 
parameters using the same folding procedure. Taking 
as a reference the solid black line (/? = 0.15), we can 
see that the main characteristics observed, such as the 
location and strength of the GT resonance are qualita- 
tively reproduced. In Table |n] we can see a compari- 
son between measured and calculated GT strengths with 
(3 = 0.15 in different energy regions. In the low energy 
region up to 5 MeV, we have two different evaluations of 
the GT strength measured in Ref. 1141 that compare well 
with the old measurements of Ref. [13[ • Our calculations 
overestimate clearly these measurements in this range of 
energy. Nevertheless, if we consider the range of energy 
up to 7 MeV, then the calculations produce comparable 
results. This is also true when we compare the strength 
contained between 7 and 10 MeV, where another peak is 
observed, as well as when we compare the strength con- 
tained between 10 and 20 MeV, where the GT resonance 
is located. The total strength in the whole energy range 



is also well described. Calculations with other deforma- 
tions produce peaks that are displaced, but qualitatively 
they are similar with one broad peak centered at about 
5 MeV and another one centered beyond 10 MeV. 

Similarly, we can see in Fig.[5]the correspondin g re sults 
for the charge-exchange reactions 76 Se(n,p) 76 As and 
76 Ge(d, 2 He) 76 As with an improved energy resolu- 
tion of 120 keV. The curves represent again folded dis- 
tributions. The lower panel shows the QRPA results for 
different values of quadrupole deformations. We see in 
this case an increased sensitivity to deformation in the 
low energy region where a first peak carrying most of 
the strength is particularly enhanced for prolate defor- 
mations. In any case, we should keep in mind that in the 
case of B(GT + ) the total strength is much lower than in 
the case of £?(GT~), as it must be according to the Ikeda 
sum rule for N > Z nuclei. Taking in this case /3 = 0.2 
(solid black line) as a reference, we can see in Table IIIII 
the strength contained in various energy ranges compared 
to experiment. The GT strength contained below 2 MeV 
is largest for the most recent data of Ref. [l6[. However, 
as in the case of 76 Ge, for energies below 5 MeV both 
sets of data and the theoretical results are comparable. 
The total strength contained in the whole energy range 
is also well reproduced by the calculations. 

This comparison tells us that the global behavior of the 
GT strength distribution as a whole is well reproduced 
in our calculations for reasonable values of the coupling 
constants of the residual interaction and for deformation 
parameters compatible with experiment. 

In the next two figures, Figs. [3] and 0] we can see the 
comparison between the QRPA results for various de- 
formations (blue lines) and the high resolution measure- 
ments for 76 Ge and 76 Se, respectively. In both figures 
we show the folded measured (solid lines) and calculated 
(dashed lines) distributions for a better comparison. 

In the case of 76 Ge we can see that we get systemati- 
cally less (more) strength than experimentally observed 
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FIG. 1: (Color online) Gamow-Teller fl(GT") strength dis- 
tributions in 76 Ge as a function of the excitation energ y in the 
daughter nucleus. The upper panel shows the data [l3l . H3 | 
from different experiments. The lower panel shows calculated 
QRPA results with various quadrupole deformations. 



below (above) an excitation energy of about 2 MeV. This 
is true for any deformation and the total GT strength 
in this region is somewhat compensated. On the other 
hand, it is clear that the highly fragmented strength mea- 
sured is only qualitatively reproduced when prolate de- 
formations of about ft ~ 0.15 are considered. In partic- 
ular, the spherical strength distribution shows a concen- 
tration of the strength in a few peaks at variance with 
experiment. The strength distributions in Fig. U corre- 
sponding to 76 Se show also the characteristic fragmen- 
tation due to deformation that agrees better with the 
experiment. In this case the total calculated strength be- 
low this small range of energy is lower than experiment 
except for large deformations that accumulate strength 
around 2 MeV (see Table EH). 

In Figs. [3] and [5] we show the accumulated GT strength 
[J2B(GT)(E ex ) = T.e<e„B(GT)(E)] for 76 Ge and 
76 Se, respectively. They are calculated from the folded 
distributions and correspond to QRPA calculations with 
various deformations using two different Skyrme forces 
to appreciate the sensitivity of the results to differ- 
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FIG. 2: (Color online) Same as in Fig.[TJ but for the B(GT + ) 
strength in 76 Se. Experimental data are from Refs. [la. HtSI]. 



ent parametrizations of the effective Skyrme interaction 
used. In the left panels (a) we have the results from SLy4, 
whereas in the right panels (b) we have the results from 
the Skyrme interaction SG2 [43] . The latter has been 
successfully tested against spin and isospin excitations in 
spherical [13] and deformed nuclei [H, HU ■ As we can see 
from these figures, the tendencies are very similar and the 
small discrepancies are only quantitative. This type of 
plot is very useful because it shows at the same time both 
the detailed structure of the strength distribution and 
the global behavior in terms of total strength contained 
at each excitation energy. In Fig. [5] the calculations for 
the accumulated B(GT~) strength in 76 Ge are compared 
to the data from 76 Ge( 3 He, t) 76 As [3] (circles) and with 
the data extracted over 0.5 MeV energy bins (triangles) 
that contain the influence of the tail of the GT resonance 
[l|j|. We also show for comparison the data from (p, n) 
reactions [Hj] (open squares). We can see that the mea- 
sured strength is systematically underestimated at low 
excitation energy, but beyond 4 MeV the tendency is the 
opposite. The data are best reproduced by the calcula- 
tions with a quadrupole deformation j3 = 0.15. In Fig. [B] 
we show the calculations for the _B(GT + ) strength in 76 Se 
for various deformations. They are compared with data 
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FIG. 3: (Color online) High resolution 76 Ge( 3 He,t) 76 As data Q3] compared to QRPA calculations (dashed and vertical lines) 
with various quadrupole deformations. 
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FIG. 4: (Color online) Same as in Fig.H but for 76 Se(d, 2 He) 76 As data [T 



from the 76 Se(<i, 2 He) 76 As [l6| (dots) and (n,p) reactions 
[l~5| (open squares). In this case the results from spher- 
ical and oblate shapes clearly underestimate the data. 
Prolate deformations produce much more strength and 
follow better the observed trend. 

Finally, it should also be mentioned that the strength 
of the residual proton- neutron interaction in the pp chan- 
nel (Kpjf) has been shown to play an important role to 
describe both the GT strength distributions JSES EH 
and the 2i//3/3 nuclear matrix elements [ol l9l. I 111 l30j|. To 
demonstrate the sensitivity of the GT strength distribu- 



tions to this parameter, we show in Fig. [7] the accumu- 
lated (a) B(GT~) and (b) B(GT + ) for several values of 
the coupling strength [MeV] . The results correspond 
to the Skyrme force SLy4 for various values of the defor- 
mation parameters close to the selfconsistent ones and to 
the experiment. We observe a dispersion of the results 
characterized by a larger accumulation of the strength at 
lower energies as the value of increases. A similar 
tendency is found for other values of the deformation pa- 
rameter. In the case of B(GT~), the spread in the profile 
produced by k£, t is comparable to the spread produced 
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E ex [MeV] E ex [MeV] 

FIG. 5: (Color online) Accumulated B(GT~) of 76 Ge as a 
function of the excitation energy in the daughter nucleus. The 
data [IH . [T3 ] are compared with theoretical calculations ob- 
tained with different quadrupole deformations. 
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FIG. 7: (Color online) Accumulated (a) B(GT~) and (b) 
B(GT + ) as a function of the excitation energy in the daugh- 
ter nucleus. The results correspond to the Skyrme force SLy4 
for various values of the coupling strength Kp^ [MeV]. Exper- 
imental points in (a) and (b) are as in Figs. [5] and [6] respec- 
tively. 




E ex [MeV] E ex [MeV] 

FIG. 6: (Color online) Same as in Fig. [5] but for the B(GT+) 
of 76 Se. Data are from Refs. [HGiJ. 



by the deformation (see Fig. Ufa)). On the other hand, 
the effect of k^J in the accumulated _B(GT + ) at low en- 
ergies is much smaller than the effect produced by the 
deformation (see Fig. HJa)). 




E ex [MeV] 

FIG. 8: (Color online) Experimental running sum [3| of the 
nuclear matrix element for the 2u/3j3 decay in 76 Ge as a func- 
tion of the intermediate excitation energy in 76 As, compared 
with results from calculations using different quadrupole de- 
formations for parent and daughter nuclei. The shaded area 
indicates the experimental range extracted from the experi- 
mental half-life using bare and quenched qa factors. 



B. Two- neutrino double-beta decay 

In this subsection we evaluate the 2vj3(3 matrix ele- 
ments for the decay of 76 Ge and compare them with the 
experimental information extracted from both the mea- 
sured half-life of the process [4l[ and the running sum, as 
extracted directly from the high resolution measurements 
00. 

In Fig. [5] we can see the running sum of the 2vj3(3 nu- 



clear matrix element for various combinations of parent 
and daughter deformations as a function of the excitation 
energy of the intermediate nucleus 76 As. The shaded area 
indicates the experimental range for the matrix element 
M"gt^ extracted from the experimental half-life using 
bare and quenched gA factors. Also included (red solid 
line) is the experimental running sum up to 5 MeV de- 
rived following the procedure explained in Ref. [lij using 
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FIG. 9: (Color online) Same running sums as in Fig. [SJ but 
for different values of the coupling strength Kpp [MeV]. 
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FIG. 10: Nuclear matrix element for the 2z//3/3 decay in 76 Ge 
as a function of the coupling strength Kpjf . The shaded area 
has the same meaning as in Fig. [8] 



the GT strength distributions from [3] and fl6l ]. 

The spherical description of these nuclei produce a very 
large matrix element that overestimates the experiment. 
The matrix elements are reduced when different defor- 
mations from parent and daughter are considered. In 
particular, the combination of deformation parameters 
/3( 76 Gc) = 0.15 and /3( 76 Se) = 0.2 (solid black line) pro- 
duce optimal results reproducing the experiment. When 
the deformations are very different, the matrix element 
becomes very small as compared to the experiment. 

It should also be noted that all the calculations lie 
below the experimental running sum of Ref. [l4| that 



reaches the experimental value from the half-life already 
at 2 MeV. As pointed out by the authors in that refer- 
ence, the consequence of this is that any further contribu- 
tions must be very small or must cancel each other. How- 
ever, one should also take into account that the construc- 
tion of the running sums in |14j was based on summing 
the products of GT - and GT + matrix elements accumu- 
lated in different overlapping energy windows. This was 
done because the different resolutions of the experiments 
on the two branches and the lack of correlation between 
the GT~ and GT + strengths do not allow a one-to-one 
correspondence between the states in 76 As excited from 
76 Ge and 76 Se. Then, this procedure provides an upper 
limit of the actual running sum because the product of 
accumulated strengths is always larger than the sum of 
one-to-one products. In this sense, our calculations for 
the running sums are consistently below the experimental 
running sum extracted in that way. 

The dependence of the running sum with the coupling 
strength of the proton-neutron residual interaction in the 
pp channel (n^p) is shown in Fig. [9] using the defor- 
mations /3( 76 Ge)=0.15 and /3( 76 Se)=0.20 for the parent 
and daughter nuclei, respectively. The concentration of 
the strength at low energies produced with higher values 
of makes the 2^/3/3 nuclear matrix element increase 
at low energies with increasing values of but nev- 

ertheless, this increase is not enough to reproduce the 
experimental running sum extracted in Ref. (l4j . When 
the running sum exhausts, the final matrix element de- 
creases with increasing values of k^ t . This effect can be 

better appreciated in Fig. [TU1 where M^ 13 is plotted as 
a function of . 



IV. SUMMARY AND CONCLUSIONS 

In this work we have studied the GT strength distri- 
butions in the daughter nucleus 76 As reached from both 
76 Ge and 76 Se double-/? decay partners. Calculations 
from a deformed QRPA approach with ph and pp residual 
interactions based on a selfconsistcnt Skyrme Hartree- 
Fock mean field with pairing correlations are compared 
with data from (p,n) and (n,p) charge- exchange reac- 
tions and their associated high resolution ( 3 He, t) and 
(<i, 2 He) reactions. 

Using quadrupole deformations compatible with the 
equilibrium shapes obtained from the SLy4 interaction, 
which lie within the experimental values, we obtain rea- 
sonable agreement with experiment in both single-beta 
branches, GT~ in 76 Ge and GT + in 76 Se, as measured in 
Refs. [IUdll and [THEHi respectively, as well as with the 
nuclear matrix element of the 2^/3/3 process, extracted 
from the experimental half-life [4l[ . The deformed QRPA 
approach used in this work is able to reproduce the gross 
features of the GT strength distributions, but fails to 
account for a detailed description of the high resolution 
data in the low excitation energy. Although a fine spec- 



9 



troscopy is beyond the scope of this HF+BCS+QRPA 
approach, where the use of universal effective Skyrme in- 
teractions (SLy4 in this work) intended to describe spher- 
ical and deformed nuclei all along the nuclear chart pre- 
vents reproducing accurately local details, the global per- 
formance of this approach is very reasonable. It is also 
worth mentioning that the agreement with experiment is 
optimum for the most reasonable choices of deformations 
and strengths of the residual interaction. 

It has been shown that nuclear deformation plays a 
significant role in understanding the GT strength distri- 
bution, as well as in understanding the 2vf3f3 process, 
where differences between parent and daughter nuclear 
deformations introduce a reduction factor in the nuclear 
matrix elements that finally determines the 2vf3/3 half 
lives. 

The role of the coupling strength of the proton- neutron 
residual interaction in the pp channel to describe 



GT strength distributions and 2vf3f3 nuclear matrix ele- 
ments has been studied. It has been shown that, within 
a range of reasonable values of k^, t , its effect on the 
GT strength distributions at low energy is comparable 
or smaller than the effect produced by the deformation. 
It has also an effect on the 2v(3/3 matrix elements, espe- 
cially on the total matrix element that finally determines 
the 2v(3j3 half-life. 
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